Error Estimates for Generalized Barycentric Interpolation 

Andrew GillettjQ, Alexander Rane0, Chandrajit BajaJ] 
April 19, 2011 

Abstract We prove the optimal convergence estimate for first order interpolants used 
in finite element methods based on three major approaches for generalizing barycentric 
interpolation functions to convex planar polygonal domains. The Wachspress approach 
explicitly constructs rational functions, the Sibson approach uses Voronoi diagrams 
on the vertices of the polygon to define the functions, and the Harmonic approach 
defines the functions as the solution of a PDE. We show that given certain conditions 
on the geometry of the polygon, each of these constructions can obtain the optimal 
convergence estimate. In particular, we show that the well-known maximum interior 
angle condition required for interpolants over triangles is still required for Wachspress 
functions but not for Sibson functions. 
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1 Introduction 

While a rich theory of finite element error estimation exists for meshes made of tri- 
angular or quadrilateral elements, relatively little attention has been paid to meshes 
constructed from arbitrary polygonal elements. Many quality-controlled domain mesh- 
ing schemes could be simplified if polygonal elements were permitted for dealing with 
problematic areas of a mesh and finite element methods have been applied to such 
meshes [351139] . Moreover, the theory of Discrete Exterior Calculus has identified the 
need for and potential usefulness of finite element methods using interpolation methods 
over polygonal domain meshes (e.g. Voronoi meshes associated to a Delaunay domain 
mesh) |18] . Therefore, we seek to develop error estimates for functions interpolated 
from data at the vertices of a polygon Q. 

Techniques for interpolation over polygons focus on generalizing barycentric coor- 
dinates to arbitrary n-gons; this keeps the degrees of freedom associated to the vertices 
of the polygon which is exploited in nodal finite element methods. The seminal work of 
Wachspress [37] explored this exact idea and has since spawned a field of research on 
rational finite element bases over polygons. Many alternatives to these 'Wachspress co- 
ordinates' have been defined as well, including the Harmonic and Sibson interpolants. 
To our knowledge, however, no careful analysis has been made as to which, if any, 
of these interpolation functions provide the correct error estimates required for finite 
element schemes. 

We consider first-order interpolation operators from some generalization of barycen- 
tric coordinates to arbitrary convex polygons. A set of barycentric coordinates {A^} for 
Q associated with the interpolation operator I : H 2 ((l) — > spanjA^} C H^{fi) is given 
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by 

Iu := y^M(vj)A». (1) 

i 

Since barycentric coordinates are unique on triangles (described in Section [5TTJ this is 
merely the standard linear Lagrange interpolation operator when fl is a triangle. 

Before stating any error estimates, we fix some notation. For multi-index a = 
(01,02) and point x = (x,y), define x Q := x ai y a ' 2 , a\ := 0102, M -~ a l + a 2, an d 
D"m := d^ a ^u/dx ai dy" 2 . The Sobolev semi-norms and norms over an open set Q are 
defined by 



\ u \h™(Q)-= I Yl l° a?t ( x )| 2dx and H|jrm (12 ) := Yl ■ 

|a|=m 0<fc<m 

The ff°-norm is the L 2 -norm and will be denoted IMI^a^y 

Analysis of the finite element method often yields bounds on the solution error in 
terms of the best possible approximation in the finite-dimensional solution space. Thus 
the challenge of bounding the solution error is reduced to a problem of finding a good 
interpolant. In many cases Lagrange interpolation can provide a suitable estimate which 
is asymptotically optimal. For first-order interpolants that we consider, this optimal 
convergence estimate has the form 

Ib'--Nlffi(fi) <Cdia,m(n)\u\ H2{a) , VuGiJ 2 (fi). (2) 

To prove estimate (0 in our setting, it is sufficient (see Section Q to restrict the 
analysis to a class of domains with diameter one and show that J is a bounded operator 
from H 2 (tt) into H 1 ^), that is 

P«lljri(n) <Ci\H\ H 2 in) , VueH\n). (3) 

We call equation Q the .H^-interpolant estimate associated to the barycentric 
coordinates A,; used to define /. 

The optimal convergence estimate © does not hold uniformly over all possible 
domains; a suitable geometric restriction must be selected to produce a uniform bound. 
Even in the simplest case (Lagrange interpolation on triangles), there is a gap between 
geometric criteria which are simple to analyze (e.g. the minimum angle condition) and 
those that encompass the largest possible set of domains (e.g. the maximum angle 
condition) . 

This paper is devoted to finding geometric criteria under which the optimal con- 
vergence estimate ([2]) holds for several types of generalized barycentric coordinates on 
arbitrary convex polygons. We begin by establishing some notation (shown in Figure 
[1]) to describe the specific geometric criteria. 

Let J? be a convex polygon with n vertices. Denote the vertices of Q by and the 
interior angle at Vj by Pi. The largest distance between two points in Q (the diameter 
of Q) is denoted diam(J7) and the radius of the largest inscribed circle is denoted p(O). 
The center of this circle is denoted c and is selected arbitrarily when no unique circle 
exists. The aspect ratio (or chunkiness parameter) 7 is the ratio of the diameter to 
the radius of the largest inscribed circle, i.e. 

diam(fi) 

7:= ' 



Fig. 1 Notation used throughout paper. 



We will consider domains satisfying one or more of the following geometric condi- 
tions. 

Gl. Bounded aspect ratio: There exists 7* G R such that 7 < 7*. 
G2. Minimum edge length: There exists d* € R such that |vj —~Vj\ > d* > for 
all i / j. 

G3. Maximum interior angle: There exists /3* £ R such that fii < f3* < n for all i. 

Using several definitions of generalized bary centric functions from the literature, 
we show which geometric constraints on fl are either necessary or sufficient to ensure 
the estimate for each definition. The main results of this paper are summarized by the 
following theorem and Table \T\ Primary attention is called to the difference between 
Wachspress and Sibson coordinates: while G[3]is a necessary requirement for Wachspress 
coordinates, it is demonstrated to be unnecessary for the Sibson coordinates. 

Theorem 1 In Table\l\ any necessary geometric criteria to achieve the H 1 mterpolant 
estimate £3|) are denoted by N. The set of geometric criteria denoted by S in each row 
are sufficient to guarantee the H 1 interpolant estimate 
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Table 1 'N' indicates a necessary geometric criterion for achieving the H 1 interpolant estimate 
J3l . The set of criteria denoted 'S' in each row, taken together, are sufficient to ensure the H 1 
interpolant estimate l|3|). 



In Section[21 we define the various types of generalized barycentric coordinates, compare 
their properties, and mention prior applications. In Section [3] we review some general 
geometric results needed for subsequent proofs. In Section [4] we give the relevant 
background on interpolation theory for Sobolev spaces and state some classical results 
used to motivate our approach. In Section 15.11 we show that the simplest technique 
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of triangulating the polygon achieves the estimate if and only if Cj3] holds. In Section 
15.21 we show that if harmonic coordinates are used, G[T] alone is sufficient. In Section 
15.31 we show that Wachspress coordinates require G[3]to achieve the estimate but all 

three criteria are sufficient. In Section 15.41 we show that Sibson coordinates achieve 

the estimate with GJT] and G[2] alone. We discuss the implications of these results and 

future directions in Section [6] 

2 Generalized Barycentric Coordinate Types 

Barycentric coordinates on general polygons are any set of functions satisfying certain 
key properties of the regular barycentric functions for triangles. 

Definition 1 Functions A, : O — > R, i = 1, . . . , n are barycentric coordinates on 

Q if they satisfy two properties. 

Bl. Non-negative: A; > on Q. 

n 

B2. Linear Completeness: For any linear function L : fi — > K, L = L(vj)Aj. 

i=i 

Remark 1 Property B2]is the key requirement needed for our interpolation estimates. 
It ensures that the interpolation operation preserves linear functions, i.e. IL — L. 

We will restrict our attention to barycentric coordinates satisfying the following invari- 
ance property. Let T : IR 2 — > M. 2 be a composition of rotation, translation, and uniform 
scaling transformations and let {\f} denote a set of barycentric coordinates on TO. 

B3. Invariance: Aj(x) = \f(T(x)). 

This assumption will allow estimates over the class of convex sets with diameter 
one to be immediately extended to generic sizes since translation, rotation and uniform 
scaling operations can be easily passed through Sobolev norms (see Section [4]). At the 
expense of requiring uniform bounds over a class of diameter-one domains rather than 
a single reference element, complications associated with handling non-affine mappings 
between reference and physical elements are avoided [3]- 

A set of barycentric coordinates {Aj} also satisfies these additional familiar prop- 
erties: 

n 

B4. Partition of unity: Aj = 1. 

i=l 
n 

B5. Linear precision: ^^VjAj(x) = x. 

i=l 

B6. Interpolation: Aj(vj) = Sij. 

The precise relationship between these properties and those defining the barycentric 
coordinates is given in the following proposition. 

Proposition 1 The properties iJ7]-iJ§| are related as follows: 
(i) 41<^ (^and 
(h) (%H and 411) => M 
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Proof Given Ej2] setting L = 1 implies E(4]and setting -L(x) = x yields E[5] Conversely, 
assuming E(4]and E(5] let L(x, y) = ax + by + c where a, 6, c £ R are constants. Let Vj 
have coordinates (vf,v^). Then 

n n 

^ L( Vi )A 4 (x-, y) = ^(avf + 6vf + c)A;(x) 

i=l i=l 

/ n \ / n \ / n 

X)vfAi(x) +& ^vfA^x) +c £>(x] 

U=l / \i=l / \i=l 

— ax + by + c — L(x, y). 

A proof that rJQand pimply rJJ]can be found in PU Corollary 2.2]. □ 

Thus, while other definitions of barycentric coordinates appear in the literature, re- 
quiring only properties E[T] and E[2] is a minimal definition still achieving all the desired 
properties. 

In the following subsections, we define common barycentric coordinate functions 
from the literature. Additional comparisons of barycentric functions can be found in 
the survey papers of Cueto et al. [H] and Sukumar and Tabarraei [M] . 



2.1 Triangulation Coordinates 

The simplest method for constructing barycentric coordinates on a polygon is to tri- 
angulate the polygon and use the standard barycentric coordinate functions of these 
triangles. Interpolation properties of this scheme are well known from the standard 
analysis of the finite element method over triangular meshes, but this construction 
serves as an important point of comparison with the alternative barycentric coordi- 
nates discussed later. 

Let T be a triangulation of Q formed by adding edges between the Wj in some 
fashion. Define 

Xjj- : n ->• K 

to be the barycentric function associated to Vj on triangles in T containing Vj and 
identically otherwise. Trivially, these functions define a set of barycentric coordinates 
on Q. 

Two particular triangulations are of interest. For a fixed i, let 7m denote any 
triangulation with an edge between v^_i and Vj + i . Let 7m denote the triangulation 
formed by connecting to all the other v^. Examples are shown in Figure [2] 

Proposition 2 (Floater et al. 110/ ) Any barycentric coordinate function Xi according 
to Definition^ satisfies the bounds 

o < x™ m (x) < \i(x) < \l r i M (x) < i, vie n. (4) 

Proposition [5] tells us that the triangulation coordinates are, in some sense, the ex- 
tremal definitions of generalized barycentric coordinates. In any triangulation of ft, at 
least one triangle will be of the form (v^_i , v^, "Vi+i), and hence the lower bound in © 
is always realized by some Xf 11 . Thus, the examination of alternative barycentric coor- 
dinates can be motivated as an attempt to find non-extremal generalized barycentric 
coordinates. 
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Fig. 2 Triangulations T m and 7m are used to produce the minimum and maximum baryccn- 
tric functions associated with Vj, respectively. 

2.2 Harmonic Coordinates 

A particularly well-behaved set of barycentric coordinates, harmonic coordinates, can 
be defined as the solutions to certain boundary value problems. Let gi : dft — > K be 
the piecewise linear function satisfying 



The harmonic coordinate function A^ ar is defined to be the solution of Laplace's equa- 
tions with gi as boundary data, 



Existence and uniqueness of the solution are well known results [13U30] . Properties E}T] 
and E[2]are a consequence of the maximum principle and linearity of Laplace's equation. 

These coordinates are optimal in the sense that they minimize the norm of the 
gradient over all functions satisfying the boundary conditions, 



This natural construction extends nicely to polytopes, as well as to a similar defini- 
tion for barycentric-like (Whitney) vector elements on polygons. Christensen [7] has 
explored theoretical results along these lines. Numerical approximations of the Ap ar 
functions have been used to solve Maxwell's equations over polyhedral grids [12] and 
for finite element simulations for computer graphics [25U21] , While it may seem ex- 
cessive to solve a PDE just to derive the basis functions for a larger PDE solver, the 
relatively limited geometric requirements required for their use (see Section [5]2} make 
these functions a useful reference point for comparison with simpler constructions and 
a suitable choice in contexts where mesh element quality is hard to control. 

2.3 Wachspress Coordinates 



9i(vj) = 5, 



gi linear on each edge of ft. 




= 0, on ft, 
= gi. on dft. 



(5) 




One of the earliest generalizations of barycentric coordinates was provided by Wach- 
spress [33 . Definition of these coordinates is defined based on some notation shown 
in Figure [3] Let x denote an interior point of ft and let Ai (x) denote the area of the 



Fig. 3 Left: Notation for A;(x). Right: Notation for Bi. 



triangle with vertices x, v.;, and v,;+i where, by convention, vo := v n and v n +i := vi. 
Let denote the area of the triangle with vertices Vj_x, v,;, and v.j + i. Define the 
Wachspress weight function 



^ Wach (x) = B t J] A,-(x). 



The Wachspress coordinates are then given by 



„,,Wacli/„\ 

,Wach r s _ ™j ( x ) / fi N 
*i W — v^n Wach/ \ \ ' 



These coordinates have received extensive attention in the literature since they can be 
represented as rational functions in Cartesian coordinates. Their use in finite element 
schemes has been numerically tested in specific application contexts but to our knowl- 
edge has not been evaluated in the general Sobolev error estimate context considered 
here. We note that A^ 11 G H 1 ^) since it is a rational function with strictly positive 
denominator on Q. 



Remark 2 Since Bi does not depend on x and ^i(x) is linear in x, the Wachspress 
functions are degree n — 2. By a result from Warren [38], the Wachspress functions 
are the unique, lowest degree rational barycentric functions over polygons. For finite 
element applications, however, the Aj need not be rational. 



2.4 Sibson (Natural Neighbor) Coordinates 

The Sibson coordinates [32, also called the natural neighbor or natural element coor- 
dinates, make use of Voronoi diagrams on the vertices v.j of O. Let x be a point inside 
Q. Let P denote the set of vertices {v, } and define 

P' = PU{x} = {vi,...,v„,x}. 

We denote the Voronoi cell associated to a point p in a pointset Q by 

V Q (p) := {y 6 n : |y - p| < |y - q| , Vq 6 Q \ {p}} . 



Fig. 4 Geometric calculation of a Sibson coordinate, d is the area of the Voronoi region 
associated to vertex Vj inside Q. D(x) is the area of the Voronoi region associated to x if it is 
added to the vertex list. The quantity -D(x) n Ci is exactly -D(x) if x = and decays to zero 
as x moves away from Vj , with value identically zero at all vertices besides v; . 



Note that these Voronoi cells have been restricted to Q and are thus always of finite 
size. We fix the notation 

Ci := |V P (vi)| = |{y G Q : |y- Vi| < |y - v, , Vj + i}\ 

= area of cell for Vj in Voronoi diagram on the points of P, 

D(x) := \V P ,{x)\ = |{y G S2 : |y - x| < |y-v. t | , Vt}| 

= area of cell for x in Voronoi diagram on the points of P' . 

By a slight abuse of notation, we also define 

£>(x)nc< - |Vp/(x)nv P (vi)|. 

The notation is shown in Figure 3] The Sibson coordinates are defined to be 

sibs, s fl(x)nCj . sib S , v D(x)nCi 

Xi (x) := — — or, equivalently, \ (x) = 



It has been shown that the Af lbs are C°° on 12 except at the vertices Vj where they 
are C° and on circumcircles of Delaunay triangles where they are C 1 [32II14| . Since the 
finite set of vertices are the only points at which the function is not C 1 , we conclude 
that Af bs G H\n). 

To close this section, we compare the intra-element smoothness properties of the 
coordinate types on the interior of S2. The triangulation coordinates are C , the Sibson 
coordinates are C , and the Wachspress functions and the harmonic coordinates are 
both C°°. 



3 Generalized Shape Regularity Conditions 

The invariance property E(3] allows estimates on diameter-one polygons to be scaled to 
polygons of arbitrary size. Several well-known properties of planar convex sets to be 
used throughout the analysis are given in Proposition [3] Let \S2\ denote the area of 
convex polygon S2 and let \dS2\ denote the perimeter of S2. 

Proposition 3 If Q is a convex polygon with diam(J?) = 1, then 
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(i) \n\ < tv /4, 

(u) \an\ < n, 

(Hi) Q is contained in a ball of radius no larger than l/y/2, and 

(iv) If convex polygon T is contained in Q, then \dT\ < \df2\. 

The first three statements are the isodiametric inequality, a corollary to Barbier's 
theorem, and Jung's theorem, respectively. The last statement is a technical result 
along the same lines. See |1 11140 lEJT] for more details. 

Certain combinations of the geometric restrictions (G1-G3) imply additional useful 
properties for the analysis. These resulting conditions are listed below. 

G4. Minimum interior angle: There exists (3* £ R such that ft > ft, > for all i. 
G5. Maximum vertex count: There exists n*£l such that n < n* . 

For triangles, G(4] and C{3] are the only two important geometric restrictions since 
G[S] holds trivially and GQj4>GU]=>G2] For general polygons, the relationships between 
these conditions are more complicated; for example, a polygon satisfying C(T]may have 
vertices which are arbitrarily close to each other and thus might not satisfy G(5] Propo- 
sition U below specifies when the original geometric assumptions (G1-G3) imply G2]or 

cm 

Proposition 4 The following implications hold. 

(i) 03^ q£ 

(ii) (CM or (J3p => CH 

Proof CQ] =>■ CQ] If ft is an interior angle, then p{Q) < sin(ft/2) (see Figure [S]). 
Thus 7 > . ,\ )„, . We conclude that ft > 2 arcsin -4- . Note that 7* > 2 so this is 

' sin(/3i/2) ' 1 7 ' — 

well-defined. 




Fig. 5 Proof that Gl => G4. The upper angle in the triangle is < ft/2 < n/2 and the 
hypoteneuse is < diam(i?) = 1. Thus p(fi) < sin(/3j/2). 



Q5]=> C0 By Jung's theorem (Proposition loTliIij) ), there exists x £ Q such that Q C 
B(x, l/y/2). By GH {B(vi,d*/2)}™ =1 is a set of disjoint balls. Thus B(x, l/y/2 + d*/2) 
contains all of these balls. Comparing the areas 

of U"=i B(vi,d*/2) and B(x, 1^/2 + 
d./2) gives < n(± + d*/2) 2 , so n < (V5 g* )2 . 

GO ^ GO Since 12 is convex, £™ =1 Pi = ~ 2 )- So ^ ~ 2 )- Thus 
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4 Interpolation in Sobolev Spaces 

Interpolation error estimates are typically derived from the Bramble-Hilbert lemma 
which says that Sobolev functions over a certain domain or class of domains can be 
approximated well by polynomials. The original lemma [5] applied to a fixed domain 
(typically the "reference" element) and did not indicate how the estimate was im- 
pacted by domain geometry. Later, a constructive proof based on the averaged Taylor 
polynomial gave a uniform estimate under the geometric restriction Cjl] |10| [6] . Recent 
improvements to this construction have demonstrated that even the condition Gjl] is 
unnecessary [361 (9] . This modern version of the Bramble-Hilbert lemma is stated below 
and has been specialized to our setting, namely, the H 1 estimate for diameter 1, convex 
domains. 

Lemma 1 ( [36j [9|) Let Q be a convex polygon with diameter 1. For all u £ H 2 {Q), 
there exists a first order polynomial p u such that \\u — Pu\\ui-(n) < Cbh \ u \h 2 (Q)- 

We emphasize that the constant Cbh is uniform over all convex sets of diameter 
1. The ii^-interpolant estimate © and Lemma [T] together ensure the desired optimal 
convergence estimate ([2]). 

Theorem 2 Let Q be a convex polygon with diameter 1. Lfthe H 1 -interpolant estimate 
0) holds, then for all u G H 2 (Q), 

\\u-lu\\ H i(n) < ( 1 + C l)\J 1 + C BH \ u \h*((2)- 

Proof Let p u be the polynomial given in Lemma [1] which closely approximates u. By 
property B^l Lp u = Pu yielding the estimate 

\\ u ~ Iu \\hi{Q) ^ \\ U -Pu\\m(n) + \\ I ( u -Pu)\\ H i {n) 

< (1 + Cr)ll«-Pu||fla(i3) < 0- + Ci)y/l + C% H \u\ H 2 (n) .□ 

Corollary 1 Let diam(i?) < 1. If the H 1 -interpolant estimate {3]) holds, then for all 
u £ H 2 (Q), 

\\ u - Iu \\HHn) ^ ( 1 + c i) \Jt + C 2 BH Ai&m(f2) \u\ H2{n] . 

Proof This follows from the standard scaling properties of Sobolev norms since property 
E(3]allows for a change of variables to a unit diameter domain. Note: the L 2 -component 
of the .H^-norm satisfies a stronger estimate containing an extra power of diam(fi). □ 

Section [5] is an investigation of the geometric conditions under which the H 1 - 
interpolant estimate © holds for the barycentric functions discussed in Section 
Under the geometric restrictions Gjl] and GO one method for verifying ([3]) (utilized 
in [6] for simplicial interpolation) is to bound the i^-norm of the basis functions. In 
several cases we will utilize this criteria which is justified by the following lemma. 

Lemma 2 Under CO] an d (3 the H 1 -interpolant estimate {3]) holds whenever there 
exists a constant C\ such that 



(7) 
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Proof This follows almost immediately from the Sobolev embedding theorem; see [2] 
[23]: 

n 
i=l 

where C s is the Sobolev embedding; i.e., |MI(7 (77) — ^ s \\ u \\H 2 (n) f° r au u e H 2 {^)- 
The constant C s is independent of the domain Q since the boundaries of all polygons 
satisfying G[T]are uniformly Lipschitz [23] . □ 



5 Error Estimate Requirements 

5.1 Estimate Requirements for Triangulation Coordinates 

Interpolation error estimates on triangles are well understood: the optimal convergence 
estimate © holds as long as the triangle satisfies a maximum angle condition [Hi 19]. In 
fact, it has been shown that the triangle circumradius controls the error independent 
of any other geometric criteria [22]. This result can be directly applied to 7 Trl , the 
interpolation operator associated to coordinates Xj* 1 . This convention will also be used 
to define I pt , J Wach y and J Slbs as the interpolation operators associated with with 
harmonic, Wachspress, and Sibson coordinates, respectively. 

Lemma 3 Under CJ3 the H 1 interpolant estimate holds for J 1 * 1 . Conversely, CJ2 
is a necessary assumption to achieve with J Tri . 

Proof All angles of all triangles of any triangulation T of Q satisfying G[3]are less than 
/?*. Thus, the sufficiency of Cj3]follows immediately from the maximum angle condition 
on simplices [3]. An example from the same paper involving the interpolation of a 
quadratic function over a triangle also establishes the necessity of the condition. □ 



5.2 Estimate Requirements for Harmonic Coordinates 

Recalling the notation from Figure [T] let T be the triangulation of Q formed by con- 
necting c to each of the ; see Figure [6^.. 

Proposition 5 Under (JJ\all angles of all triangles ofT are less than n — arcsin(l/7*). 

Proof Consider the triangle with vertices c, Vj and Vj+i. Without loss of generality, 
assume that jc — v.;| < |c — v;+i|. First wc bound Zcv.;+iVi. By the law of sines, 

sin(Zcv. t+ iv I ) = |c-Vjj < 1 
sin(Zcv l v. 1+ i) |c-v l+ i| 

If ZcVjVi_|_i > 7r/2 then ZcVj+iV; < 7r/2. Otherwise, (jSJ implies Zcv^+i Vj < tt/2. 

To bound angle ZcVjVi_|_i, it suffices to consider the case when ZcVjVj + i > n/2, 
as shown in Figure [6p. Define y to be the point on the line through v.; and v.;_|_i which 
forms a right triangle with Vj and c. Since ZcvjVi+i > it/2, y is exterior to fi, as 
shown. Observe that 

Ic — v„-l |c — v.i i diam(fi) * 
< —, < 7~7T\ — = 7 < 7 . 



|c-y| |c-y| p{0) 
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(a) 



(b) 



(c) 



Fig. 6 (a) Triangulation used in the analysis of Harmonic coordinates, (b) Notation for proof 
of the bound for ZcviV;_|_i in a case where it is > n/2. (c) Notation for proof of the bound 
for Zvicv^i in a case where it is > n/2. 



Since sin(7r - ZcVjV l+ i) = 

For the final case, it suffices to assume ZvjCVj+i > n/2, as shown in Figure [B};. 
Define y in the same way, but note that in this case y is between Vj and v,-_|_i, as 
shown. Similarly, ^r^z^p < 7*, implying ZvjVj + ic > arcsm(l/7*). Since ZvjCVj + i < 

□ 



the result follows. 



n — Zv^Vi+ic, the result follows 



Lemma 4 Under CO] the operator J° pt satisfies the H 1 interpolant estimate fgp . 

Proof Since the differential equation (J5]) is linear, I° pt u is the solution to the differen- 
tial equation, 

= 0, on J? 
L pt u = g u , on dfl 



(9) 



where g u is the piecewise linear function which equals u at the vertices of fl. Follow- 
ing the standard approach for handling nonhomogeneous boundary data we divide: 
I pt u = Whom + Wnon where Unon £ H 1 (fl) is some function satisfying the boundary 
condition (i.e., u n0 n = gu on dfl) and Uhom solves, 



4u hom = — ZiUnon, On fl 

u hom = 0, on dfl. 



(10) 



Specifically we select u n on to be the standard Lagrange interpolant of u over triangu- 
lation T (described earlier). Since Proposition [5] guarantees that no large angles exist 
in the triangulation, the standard interpolation error estimate holds, 



Unon||jyi(fi) < C B A \ u \ll -(n) 



(11) 



where Cba only depends upon the aspect ratio bound 7* since diam(fi) = 1. The 
triangle inequality then implies that ||wnon||jji(^) < max(l, Cba) I M \h 2 (C2) ■ 

Next a common energy estimate (see [13]) for (|10p implies that |uhoml//i (o) < 



|Wnon|jji(f2) ■ The Poincare inequality (see [23]) ensures that ||u n , 



om\\ L 2 (f2) 



< C P \u\ 



ff!(J7) 
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where Cp only depends on the diameter of Q which we have fixed to be 1. The argument 
is completed by combining the previous estimates: 



rOpt 



HliPffi) - H u liomllffi(j7) + ||«non||Hi(fi) 

< (1 + C P ) | 

u lwm\ H l (S}) + 1 1 "non I \h 1 (S2) 

< (1 + Cp) \\u non \\ H1{s - 2) < (l + C P )msx(l,C BA ) \\u\\ H2{n) .□ 



5.3 Estimate Requirements for Wachspress Coordinates 



VI VI VI 




Fig. 7 Example showing the necessity of condition G(3]for attaining the optimal convergence 
estimate ^ with the Wachspress coordinates. As the shape approaches a square, the level sets 
of A^ 11 collect at the top edge, causing a steep gradient and thus preventing a bound on the 
H 1 norm of the error. The figures from left to right correspond to e values of 0.1, 0.05, and 
0.025. 



Unlike the harmonic coordinate functions, the Wachspress coordinates can produce 
unsatisfactory interpolants unless additional geometric conditions are imposed. We 
present a simple counterexample (observed qualitatively in [T5] and in Figure [7} to 
show what can go wrong. 

Let Q € be the pentagon defined by the vertices 

Vi = (0,1 + 6), V 2 = (l,l), v 3 = (l,-l), V 4 = (-l,-l), V 5 = (-l,l), 

with e > 0. As e — ► 0, O e approaches a square so Gl is not violated. Consider the 
interpolant of u(x) = 1 — x\ where x = (xi, xz)- Observe that u has value 1 at vi and 
value at the other vertices of Q e . Hence 



r Wach 
1 u ■ 



^u(v,;)A 



Wach 



,Wach 



i=l 



Using the fact that du/dy = 0, we write 



r Wach 
1 u 



2 



, Wach 



2 

\h 1 (C2. 



1 x Wach 1 2 

|u-Ai | 



d(u 



Af ach ) 



d.v 



2 


aA Wach 


+ 


dy 
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The last term in this sum blows up as e — > 0. 



Lemma 5 lim 

e->0 





g A Wach 


j 


dy 



2 

= CO. 



The proof of the lemma is given in the appendix. As a corollary, we observe that 
j | ti — Iu\\ H i^ n ^ cannot be bounded independently of e. Since IMI/p^) i s finite, this 
means the optimal convergence estimate ([2]) cannot hold without additional geometric 
criteria on the domain Q. This establishes the necessity of a maximum interior angle 
bound on vertices if Wachspress coordinates are used. 

Under the three geometric restrictions GJT] G(2j and G0 ([2]) does hold which will 
be shown in Lemma [5] We begin with some preliminary estimates. 

Proposition 6 For all x G ft, \X7Ai(x)\ < ^. 

Proof In |17l Equation (17)] it is shown that the |VAi(x)| = ^ | v i — v i+l|- Since 
diam(J2) = 1 the result follows. □ 

Next we show that the triangular areas are uniformly bounded from below given 
our geometric assumptions. 

Proposition 7 Under CJ7J and (£3 there exists B* such that Bi > B*. 

Proof By G[2] the area of the isosceles triangle with equal sides of length d* meeting 
with angle /?j at v.j is a lower bound for Bi, as shown in Figure [5] (left). More precisely, 
Bi > (d*) 2 sin03j/2) cosCSi/2). (£2 implies that cos(ft/2) > cos(/3*/2). Gg] (which 
follows from G[T] by Proposition [4| implies that sin(/?i/2) > sin(/3*/2). Thus Bi > 
B* := (d*) 2 sin(/3*/ 2 ) cos(^*/2). □ 

Proposition [7] can be extended to guarantee a uniform lower bound on the sum of 
the Wachspress weight functions. 

Proposition 8 Under Gj] (J3 an d £0 there exists ui* such that for all x £ Q, 

w k [x)>w*. 

k 

Proof Let Vj be the nearest vertex to x, breaking any tie arbitrarily. We will produce 
a lower bound on u>i(x). Let j ^ {i — G[2] implies that |x — vy > d*/2 and 

|x — Vj_)-i| > d*/2. G[3] implies that ZxvjVj+i < j3* and Zxvj+iVj < /3*. It follows 
that Aj(x.) > (d*) 2 sin(7r — /3*)/4 (see Figure [8] (right)). We now use Proposition and 
property G(S] (which follows from either or GO by Proposition g]) to conclude that 



5>f aCh (x) > Wj Wach (x) J] A,-(x) > B* [(d*) 2 sin(7r - 0*)/4] 



-2 



Lemma 6 Under Gl, G2, and G3, ^ holds for the Wachspress coordinates. 



15 



Vj_l 




Vi+l 




V J+1 



Fig. 8 Left: Justification of claim that Bi > (a!*) 2 sin(/3i/2) cos(/3i/2) in the proof of Propo- 
sition [7] The shaded triangle is isosceles with angle ft and two side lengths equal to d* as 
indicated. Computing the area of this triangle using the dashed edge as the base yields the 
estimate. Right: Justification of claim that Aj(x) > (d*) 2 sin(7r — /3*)/4 in the proof of Propo- 
sition [8] The indicated angle is at least it — /3* by G(3]and \vj — Vj_|_i| > d*. Computing the 
area of the triangle using edge VjVj_|_i as the base yields the estimate. 



Proof The gradient of Af ach (x) can be bounded using Proposition [5] 



VAf ach (x) < 



|Vwf acn (x)| 



..Wach 



+ 



WE fc |v< ach W 

£,^ Wach (x)) 2 



V^ Wach (x) 


+ J2 k 


V W Wach (x) 


2E fe 

- < 


V«^ ach (x) 


Ej w T ach 


(x) 







(12) 



Recalling Proposition [3] Ej=i A;( x ) < 7r /4 and < 7r/4. Using Proposition [6] and 
the arithmetic mean-geometric mean inequality, we derive 



y - 

^ 8 



BilVA^x)) J] A fe (x 

Z)fc^i-l,i,j ^fc(x) 



jjti-l,i 



= -In -2) 
8 y ' 



71 — 3 
-i n— 3 



- ^ 8 



4(n - 3) 



n— 3 



4(n - 3) 



(13) 



By induction, one can show that n(n — 2) ^ 4(7^-3) j — ^ n ^ or n — 4- Using this, we 
substitute |T3J) into (TT2")) to get 



n— 3 



|vAf ach (x)| < A^| Vro W ach (x)| < A nj (n-2) 



4(n - 3) 



-I n-3 2 
< JL 2 ,T=^. 

4to* 2ui* 
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Fig. 9 Notation for proof of Proposition [9] 



Since \ fl\ < 7r/4 by Proposition we thus have a uniform bound 



,Wach 



4w? / 4 



5.4 Estimate Requirements for Sibson Coordinates 

The interpolation estimate for Sibson coordinates is computed using a very similar 
approach to that of the previous section on Wachspress coordinates. However in this 
case the geometric condition G[3] is not necessary. We begin with a technical property 
of domains satisfying conditions G[T]and G[2] 

Proposition 9 Under d[J\and G{21 there exists h* > such that for all x £ O, B(x, h*) 
does not intersect any three edges or any two non-adjacent edges of Q. 

Proof Let x £ O, h £ (0, d*/2), and suppose that two disjoint edges of Q, ej and ej, 
intersect B(x,h). Let Li and Lj be the lines containing e, and ej and let be the 
angle between these lines; see Figure [9] We first consider the case where Lj and Lj are 
not parallel and define z = L j n L j . 

Let Vj and Vj be the endpoints of and ej nearest to z. Since h < d*/2 both 
and vj cannot live in B(x, h); without loss of generality assume that v.j ^ B(x, h). 

Since dist(vj ■, Lj) < 2/i, 

sin0<2/i/|z-Vj|. (14) 

Let W be the sector between Li and Lj containing x. Now 12 C B(vj, 1) n W C 
B(z, 1 + |z - v, |) n W. It follows that p(J?) < (1 + |v, - ssl) sin0. Using dHJ and GQ 



< 



2/i 



(1+ z - VJ 



where the final inequality holds because by G[2] I z — Vj\ > Vj — v^- > d* . Thus 



h > 



2 7 *(l + d*) 



(15) 



Estimate ()15|l holds in the limiting case: when Li and L^ are parallel. In this case Q 
must be contained in a strip of width 2h which for small h violates the aspect ratio 
condition. 

The triangle is the only polygon with three or more pairwise non- adjacent edges. 
So it remains to find a suitable h* so that -B(x, h*) does not intersect all three edges of 
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Fig. 10 The proof of Proposition HOI has two cases based on whether x is within h„/2 of some 
Vi or not. When x is within h*/2 of Vi, the shaded sector shown on the right is contained in 
V P , (x) n O. 

the triangle. For a triangle, p(H) is the radius of the smallest circle touching all three 
edges. Since under GQ] p(ft) > 1/7*, B(x, ^r) intersects at most two edges. Thus 
ft* = 2-y*(i+d ) * s sufficiently small to satisfy the proposition in all cases. □ 

Proposition [5] is a useful tool for proving a lower bound on -D(x), the area of the 
Voronoi cell of x intersected with Q. 

Proposition 10 Under G[7] and there exists _D* > such that D(x) > D* . 

Proof Let ft* be the constant in Proposition [5] We consider two cases, based on whether 
the point x is near any vertex of O, as shown in Figure [10] (left). 
Case 1 : There exists Vj such that x £ B(yi, ft*/2). 

Consider the sector of B(x, ft*/2) specified by segments which are parallel to the 
edges of fi containing Vj, as shown in Figure [TO] (right). This sector must be contained 
in n by Proposition [9] and in the Voronoi cell of x by choice of ft* < rf* . Thus by Gj4] 
(using Proposition gjll)) D(x) > /3*ft*/8. 
Case 2 : For all v>, x £ S(vj, ft*/2). 

In this case, B(x, ft* /4) n Q C Vp' (x). If -B(x, ft*/4) intersects zero or one bound- 
ary edge of Q, then D(x) > 7tft*/32. Otherwise _B(x, ft*/4) intersects two adjacent 
boundary edges. By GH D(x) > /3*ft^/32. □ 

General formulas for the gradient of the area of a Voronoi cell are well-known and 
can be used to bound the gradients of D(x) and -D(x) n C;. 

Proposition 11 |VD(x)| < tt and |V(-D(as) n Q)| < 1. 

Proof The gradient of the area of a Voronoi region is known to be 

n 

, , \ ^ V ,- — X 

VC(X)=^^ rFj, 

3=1 1 J 

where Fj is the length of the segment separating the Voronoi cells of x and Vj [27 , 28 . 
Then applying Proposition [3] gives 

n 

|VD(x)| < J^Fi < \dO\ < tt. 

1=1 
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Similarly, 

V(-D(x) l~l Ci) = 7-^3— r -^i, 

and since F t < diam(fi), |V(D(x) n C l )\ < 1. □ 

Propositions [10] and [TTJ give estimates for the key terms needed in proving @ for 
the Sibson coordinates A Slbs . 

Lemma 7 Under Gl and G2, £7]) holds for the Sibson coordinates. 
Proof |\7Af lbs | is estimated by applying Propositions HOI and 1111 

UjSibel < |V(g(x)nCQl (Q(x)nc t )lVD(x)i [vp(x) n d)\ + [V£>(x)| 

I ' 1" D(x) + L>(x) 2 " D(x) 

< 1 + 7r 

Integrating this estimate completes the result. □ 

Corollary 2 By Lemma\^ the H 1 interpolant estimate holds for the Sibson coor- 
dinates. 

6 Final Remarks 

Geometric requirements needed to ensure optimal interpolation error estimates are 
necessary for guaranteeing the compatibility of polygonal meshes with generalized 
barycentric interpolation schemes in finite element methods. Moreover, the identifica- 
tion of necessary and unnecessary geometric restrictions provides a tool for comparing 
various approaches to barycentric interpolation. Specifically we have demonstrated the 
necessity of a maximum interior angle restriction for Wachspress coordinates, which 
was empirically observed in [TS], and shown that this restriction is unneeded when 
using Sibson coordinates. 

Table [TJ provides a guideline for how to choose barycentric basis functions given 
geometric criteria or, conversely, which geometric criteria should be guaranteed given 
a choice of basis functions. While utilized throughout our analysis, the aspect ratio 
requirement GQ]can likely be substantially weakened. Due to a dependence on specific 
affine transformations, such techniques on triangular domains [4l|19j (i.e., methods for 
proving error estimates under the maximum angle condition rather than the minimum 
angle condition) cannot be naturally extended to polygonal domains. Although aimed 
at a slightly different setting that we have analyzed, challenges in identifying sharp ge- 
ometric restrictions are apparent from the numerous studies on quadrilateral elements, 
e.g., [201141 |[Tll24| . A satisfactory generalization of the maximum angle condition to 
arbitrary polygons is a subject of further investigation. 

This paper emphasizes three specific barycentric coordinates (harmonic, Wachs- 
press, and Sibson) but several others have been introduced in the literature. Maximum 
entropy [33], metric [26], and discrete harmonic [29] coordinates can all be studied 
either by specific analysis or generalizing the arguments given here to wider classes of 
functions. The mean value coordinates defined by Floater [TF] are of particular interest 
in this regard as they are defined by an explicit formula and appear to not require a 
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maximum angle condition. The formal analysis of these functions, however, is not triv- 
ial. Additional generalizations could be considered by dropping certain restrictions on 
the coordinates, such as non-negativity, or the mesh elements, such as convexity. Work- 
ing with non-convex elements, however, would require some non-obvious generalization 
of the geometric restrictions G1-G5. 



A Proof of Lemma l5l 



Proof The explicit formula for the Wachspress weight associated to vi is 

w Wach (x) = BlA2 A 3 A 4 = e (l _ x )(l + x )(l + y) 

where x = (x, y) is an arbitrary point inside Q t . The other weights can be computed similarly, 
yielding the coordinate function 



I0j(v) 



e(l-x)(l + x)(l + y) 



E?=i»»j(v) e 2 (l-x 2 ) + 4e + 2(l~y) 
The y partial derivative term is computed to be 

aA Wach ^ 4e ^ _ x 2j + 4e 2(! _ x 2) + £ 3^ _ ^2)2 
dy ~ {e 2 (l-x 2 )+Ae + 2{l-y)) 2 

Define the subrcgion Q' p C fit by 

Q' p = Ux,y) e n e : i < x < ~, 1 < y < 1 + ej 
Observe that j~ < 1 — x 2 < y2 on Q' p . Fix < e < 1. We bound the numerator by 



4e (l-x 2 )+4e 2 (l-x 2 ) + e 3 (l-x 2 ) 2 > 4e ■ — + 4e 2 ■ — + e 3 ■ — > -e. 
v ' K ' K ' 16 16 256 4 

Since \y — 1| < e on fi' p , we can bound the denominator by 

|c 2 (l - x 2 ) + 4e + 2(1 - y)\ < |e 2 (l - x 2 )\ + |4e| + |2(1 - y)\ < e 2 + At + 2e < 7e. 
Putting these results together, we have that 



AL = ±_ 
49e 2 28e 



> 0. 



Let C = ^ for ease of notation. Since \Q' p \ > ie, 



lim / i > lim / 

^0j n dy ~^0j n , 



> lim 



°Jsi> 



C 2 hm^>^ 
e^O e 2 8 



lim — = co, 
e->0 e 



thereby proving the lemma. 
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